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Abstract 

The weak Stratonovich integral is defined as the limit, in law, of Stratonovich-type 
symmetric Riemann sums. We derive an explicit expression for the weak Stratonovich 
integral of f(B) with respect to g(B), where B is a fractional Brownian motion with 
Hurst parameter 1/6, and / and g are smooth functions. We use this expression to 
derive an Ito-type formula for this integral. As in the case where g is the identity, the 
Ito-type formula has a correction term which is a classical ltd integral, and which is 
related to the so-called signed cubic variation of g(B). Finally, we derive a surprising 
formula for calculating with differentials. We show that if dM = X dN, then Z dM 
can be written as ZX dN minus a stochastic correction term which is again related to 
the signed cubic variation. 

AMS subject classifications: Primary 60H05; secondary 60G15, 60G18, 60G22. 

Keywords and phrases: Stochastic integration; Stratonovich integral; fractional 
Brownian motion; weak convergence. 

1 Introduction 

If X and Y are stochastic processes, then the Stratonovich integral of X with respect to Y 
can be denned as the ucp (uniformly on compacts in probability) limit, if it exists, of the 
process 

:) • -V(/,: 



-(y^-yCVO), 



ti<t 



as the mesh of the partition {tj} goes to zero. If we specialize to the uniformly spaced 
partition, tj = j/n, then we are interested in the Stratonovich-type symmetric Riemann 
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sums, 

E x( H +xfe V fe)-yfe-0), (i.i) 

where |_^J denotes the greatest integer less than or equal to x. 

It is well-known (see [2] and [I]) that if Y — B H , a fractional Brownian motion with 
Hurst parameter H, and X = f(B H ) for a sufficiently differentiable function /, then the 
Stratonovich integral of X with respect to Y exists for all if > 1/6, but does not exist for H = 
1/6. Moreover, if H > 1/6, then the Stratonovich integral satisfies the classical Stratonovich 
change-of-variable formula, which corresponds to the usual fundamental theorem of calculus. 

In [6], we studied the case H = 1/6. There we showed that if Y = B = B 1 / 6 and 
X = f(B), where / G C°°(R), then the sequence of processes (II. ip converges in law. We let 
J*f{B(s))dB(s) denote a process with this limiting law, and we referred to this as the weak 
Stratonovich integral. We also showed that the weak Stratonovich integral with respect to 
B does not satisfy the classical Stratonovich change-of-variable formula. Rather, it satisfies 
an Ito-type formula with a correction term that is a classical It 6 integral. Namely, 

f(B(t)) = f(B(0)) + jf /'(B(a)) dB(s) - 1 jf f'"(B(s)) d[B]» (1.2) 

where [B] is what we called the signed cubic variation of B. That is, [B] is the limit in law of 
the sequence of processes ^2^l(B(tj) — B(£j_i)) 3 . It is shown in [7] that [B] = kW, where 
W is a standard Brownian motion, independent of B, and k is an explicitly defined constant 
whose approximate numerical value is k ~ 2.322. (See (12.11) for the precise definition of 
k.) The correction term above is a standard Ito integral with respect to Brownian motion. 
Similar Ito-type formulas with an Ito integral correction term were developed in [1] and [5]. 
There, the focus was on quartic variation processes and midpoint- style Riemann sums. A 
formula similar to (II. 2ft . but with an ordinary integral correction term, was established in 
[3] for the Russo-Vallois symmetric integral with respect to finite cubic variation processes. 

The precise results in [T] and [6], as well as in this paper, involve demonstrating the 
joint convergence of all of the processes involved, with the type of convergence being weak 
convergence as processes in the Skorohod space of cadlag functions. In Section [21 we establish 
the formal definition of the weak Stratonovich integral as an equivalence class of sequences 
of cadlag step functions, and we demonstrate in Theorem 12.31 the joint convergence in law of 
such sequences. For simplicity, we omit discussion of these details in this introduction, and 
only summarize the results of Section [3j in which we derive our various change-of-variable 
formulas. 

In Section [3], we extend the Ito-type formula ( 11. 2ft to the case Y = g(B). We show that 
the sequence of processes ( 11.11) converges in law to an integral satisfying the Ito-type formula 

<p(X(t)) = <p(X(0)) + jf v'(Y(s)) dY(s) - 1 jT cp"'(Y(s)) d[Y%, (1.3) 

where 

m= [\g'(B(s))fdlB} s 
Jo 
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is the limit, in law, of X/|=i(^"(^i) ~ ^(^?-i)) 3 - That is, [Y] is the signed cubic variation of 

y. 

This result is actually just one of the two main corollaries of our central result. (See 
Corollary 13.61 ) To motivate the other results, consider the following. Formulas such as 
(11. 2h and ( 11.31) are typically referred to as change-of-variable formulas. They have the same 
structure as Ito's rule, which is also generally referred to as a change-of-variable formula. 
In elementary calculus, we perform a change-of-variable when we convert an integral with 
respect to one variable into an integral with respect to another. In Ito's stochastic calculus, 
we may wish to convert an integral with respect to one semimartingale into an integral with 
respect to another. Strictly speaking, Ito's rule is not sufficient for this purpose. Ito's rule 
simply tells us how to expand a function of a semimartingale into a sum of integrals. In 
order to convert one integral into another, we must combine Ito's rule with a theorem that 
says: 

IfM = JxdY, then J Z dM = J ZXdY. 

Or, in differential form: 

If dM = X dY , then Z dM = ZX dY . (1.4) 

For It 6 integrals, this theorem is usually proved very early on in the construction of the 
integral. It is also true for the classical Stratonovich integral for semimartingales, as well 
as for ordinary Lebesgue-Stieltjes integrals. In fact, in the theory of Lebesgue-Stieltjes 
integration, it is often this result which is called the change-of-variable formula. 

In terms of the calculus of differentials, Ito's rule tells us that if M — f{Y), then 
dM = f(Y)dY + \f"{Y)d{Y), where (Y) is the quadratic variation of Y; and (Oj) 
tells us that it is permissible to substitute this expression into ZdM, so that ZdM = 
Zf{Y)dY + \Zf"{Y)d{Y). 

In this paper, we will show that (II. 4p is not true for the weak Stratonovich integral. A 
very simple example which illustrates this is the following. First, let us note that when the 
integral is defined as a limit of Stratonovich-type symmetric Riemann sums, it is always 
the case that J 6 d6 = |# 2 , for any process 9. Let us therefore define M = \B 2 , so that 
AM = B dB. On the other hand, 

M dM = -M 2 = -B 4 . 
2 8 



ifi 3 dfi-_L J W dlB] = J MBdB- 1 - J Bd[B]. 



Using (LL2J), we have 
8 

It follows that, in this example, f ll.4j) does not hold for the weak Stratonovich integral. 
Instead, we have that dM = B dB, whereas M dM = MB dB - \B d{B\. 

The second main corollary of our central result is that the weak Stratonovich integral 
satisfies a rule analogous to (II. 4p . but with a correction term. (See Corollary 13.71 ) Namely, 
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suppose X = f(B), Y = g(B), and Z = h(B), where f,g,h G C°°(R). Then the weak 
Stratonovich integral satisfies the following rule for calculating with differentials: 

If dM = XdF, then Z dM = ZX dY - -{fg'h'){B) d{B\. (1.5) 

We actually prove a slightly more general rule; see (j3.4p . 

Both p.3p and (jl.5p will be demonstrated as corollaries of the following general result. 
With X and Y as above, 

j\(s) dY{s) = *(B(t)) - $(5(0)) + 1 j\f"g' - fg")(B(s)) d[B]„ (1.6) 

where $ G C°°(R) is chosen to satisfy $' = /</. See Theorem 13.31 for the precise statement. 
Theorem 13.31 is actually formulated more generally, for integrators of the form Y + V, where 
V = J 0(B) d[-B]. This generalization is necessary to make sense of J Z dM in (ll.5p . since 
if M = J X dy, then according to (11. 6p . M is not a function of 5, but is rather the sum of 
a function of B and a process V which is in an integral against [B\ . 

2 Notation and definitions 
2.1 Basic notation 

Let B = B 1 / 6 be a fractional Brownian motion with Hurst parameter H = 1/6. That is, B 
is a centered Gaussian process, indexed by t > 0, such that 

£[£(s)£(t)] = i(t 1/3 + s 1/3 - \t - si 1 / 3 ). 

For compactness of notation, we will sometimes write B t instead of B(t), and similarly for 
other processes. Given a positive integer n, let tj = tj n = j/n. We shall frequently have 
occasion to deal with the quantity 

Let A_B jV „ = B(tj) - Bitj^) and fi*(T) = sup < t < T \B(t)\. 
Let k > be defined by 

^ = I J^dr + 1| V3 + \r - 1| 1/3 - 2|rf/ 3 ) 3 . (2.1) 

rGZ 

Let D R d[0, oo) denote the Skorohod space of cadlag functions from [0, oo) to M. d . Throughout 
the paper, "^>" will denote convergence in law. The phrase "uniformly on compacts in 
probability" will be abbreviated "ucp." If X n and Y n are cadlag processes, we shall write 
X n fa Y n or X n (t) w Y n (t) to mean that X n — Y n ucp. 
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2.2 The space [S 

Recall that for fixed n, we defined tk = k/n. Let S n denote the vector space of stochastic 
processes {L(t) : t > 0} of the form L = J2T=o ^k^[t k ,t k+1 ), where each A*,. G J 7 ^. Note that 
Xk = L{tk). Given L G S n , let Sj(L) = L{tj) — L(tj-i), for j > 1. Since t G [tk,tk+i) if and 
only if \nt\ = k, we may write 

[nt\ 

L(f) = L(0) + J>(L). 
j=i 

Definition 2.1. Let S denote the vector space of sequences A = {A n }^L 1 such that 

(i) A n G S n , 

(ii) A n (0) converges in probability, and 

(iii) there exist ipi,(p 3 ,ip 5 G C°°(M) such that 

5,{A n ) = </?iC5,-)A£,- n + tpzifoABln + (p 5 ((3j)ABj jTl + R hn , (2.2) 
where for each T, K > 0, there exists a finite constant Ct,k such that 

|-Ri,n|l{_B*(T)</<} < Ct,kt| A_Bj !rl | 7 , 

whenever j/n < T. 
If X = /(B), where / G C°°(R), then we define 

oo 

A n = JI- X '(*fc) 1 [tik,tfc +1 )> 
fc=0 

and A x = {A*}^. Note that the map X ^ A x is linear. 

Lemma 2.2. If X = f(B), where f G C°°(R) ; then A x G S and A x X uniformly on 
compacts a.s. 

Proof. Since X is continuous a.s., we have that A x — > X uniformly on compacts a.s. 
Clearly, A* G S n and A*(0) = X(0) for all n, so that Definition 0(i) and (ii) hold. For 
a, ft G R, we use the Taylor expansion 

f(b) - f(a) = f'(x)(b -a) + ±f m (z)(b - a) 3 + ^-J {5) (x)(b - af + h(a, b)(b - a) 7 , 

where x = (a + b)/2, and \h(a, b)\ < M(a,b) = sup a . 6 [ aA6)(lV6 ] \g^(x)\. For a derivation of this 
Taylor expansion, see the proof of Lemma 5.2 in [6]. 
Taking a = B(tj-i) and b = B(tj) gives 



tyAj) = f(B(tj)) - f(B(tj-i)) 

(2.3) 

24' hn 5!2 4 ' 



/'()%) AS i>B + -/"' (ft) ABJ >B + ^/ (5) (^A4 + % 



where \R j>n \ < M(B(t j - 1 ), B(tj))\AB jtn \ 7 . lij/n < T and B*(T) < K, then B(tj) G 

[-K,K], which implies M{B{t j - 1 ),B{t i )) < sup xg j_^- ^ | g^ (x) | < oo, and this verifies 
Definition I2.1( iii) showing that A x G S. □ 

We may now identify X = f(B) with A x G 5, and will sometimes abuse notation by 
writing X G S. In this way, we identify the space of smooth functions of B with a space 
of sequences in such a way that each sequence converges a.s. to its corresponding process. 
What we see next is the every sequence in S converges to a stochastic process, at least in 
law. 

Theorem 2.3. Let A (1) , . . . , A (m) G S. For 1 < k < m, choose ip lsk: <p 3jk , <p 5>k G C°°(R) 
satisfying (12. 2p for A^ and let I^ k \0) be the limit in probability of An (0) as n — >• oo. Lei 
£ C°°(IR) satisfy = and $fc(0) = 0. Let W be a Brownian motion independent of 
B, and let k > be given by ( 12. ip . Define 

X^\t) = J( fc )(0) + S fc (S(t)) + « jf (^3, fc - ^, fc ) (B(s)) dW(s), 

where this last integral is an ltd integral. Then (B,An\ ■ ■ ■ , An) =>- (BjX^, . . . ,X^ m - ) ) in 
Dgm+i [0, oo) as n -4- oo. 

Proof. By Definition I2.1[ we may write 

[nt] [nt] [nt] 

A«(t) = A«(0) • ^rl.lAh)^,,, +J2^M) AB ln + J2^M) AB l + Mt): 

i=i i=i i=i 

where i? n (t) = Ej=i J let R n (T) = sup < t < T \R n (t)\ < Y.fli \RjA- Let e > and 

choose # such that P(B*(T) > K) < e. Then 

P{R* n {T) >e)< P(B*(T) > K) + p(c t , k £ |A5 jin | 7 > e J . 

Since -B has a nontrivial 6- variation (see Theorem 2.11 in [6]), we have Yl^=i \ AB j,n\ 7 ~~ ^ 
a.s. Hence, for n sufficiently large, we have P(R n (T) > e) < 2e, which gives R n — > ucp. 

As in the proof of Theorem 2.13 in [6], we may assume without loss of generality that 
each ipi jk has compact support. By Lemma 5.1 in [6], if tp G C 1 (1R) has compact support, 
then £}=i ^(^)ABj )n ucp. Thus, 

Similarly, by Q, 

* fc (B(f)) » ^(^(Bfo)) - MBit^i))) 

3=1 

[nt] [nt] 



J2 <PlM)*B 3> + <M AB L- 
3=1 3=1 
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Therefore, 

\nt\ 

A«(f) « X (fc) (0) + $,(£(*)) + £ ^C«,0ABj B) 

where ^ = ^ - £< fc . Let K(</>, *) = E|=i Wilkin ™d J fc (t) = « £ 

By Lemma 5.2 and Theorem 2.13 in [6], we have (B, Ki(V'i), . . . , Ki(^m)) =>■ (B, Ji, . . . , J m ), 

in D R m+i[0, oo) as n -»■ oo, which implies (5, A^, . . . , A^ m) ) => (B,X^\ . . . ,X^). □ 

We now define an equivalence relation on S by A = if and only if A n — O n — > ucp. 

Lemma 2.4. I/A 6 5, i/ien i/iere exzsi unique functions fi,(f3 which satisfy (12. 2p . 7/ we 
denote these unique functions by (p\^\ and <^3,a ; f/ien A = Q if and only if 

(i) A n (0) — 6 n (0) — >■ in probability, and 

(ii) ^i,a = Vi,e an d ^3,a = ^3,6- 

Proof. Let A G S. Let {<^i, (^3, ^5} and {ipi, ^3, ^5} be two sets of functions, each of which 
satisfies QZ2). Let 1(0) be the limit in probability of A n (0) as n ->• 00. Let $, $ G C°°(R) 
satisfy $' = ipi, $' = <£>i, and $(0) = $(0) = 0. Then, by Theorem 12.31 A n converges in law 
in Dr[0, 00) to 

Z(f ) = X(0) + *{B(t)) +kJ" L 3 - Lp'^j (B(s)) dW(s) 

= 1(0) + £(B(t)) + K jf (V 3 - ^1) dW(fl). 

Hence, E[I(t) -1(0) \ F*] = $(B(t)) = $(B(t)) a.s. for all t > 0, which implies $ = $, 
and hence, = <^i. It follows that 

Jo 

Hence, £[.M(t) 2 | 7^] = /„* |(^ 3 - &)(B(s))\ 2 ds = a.s. for all t > 0, which implies 
ips = (f>3- This shows that there exist unique functions <£>i,a, <^3,A which satisfy (12.21) . 

Let A, O G S and define T = A — O. Note that A n — n — > ucp if and only if T n ^> 
in Z? R [0, 00). 

First assume (i) and (ii) hold. Then r„(0) — > in probability, so by Theorem 12.3} T n 
converges in law in Dr[0, 00) to 

MB(t)) +k£ Ls,v ~ ^<r) (B(s)) dW(s), 

where $p = <pi t r and $r(0) = 0. But from (12.21) . we see that <pip = <fi,A — fi,e = and 
V?3,r = <f3,A — y?3,e = 0. Hence, r n =>■ and A = 0. 
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Now assume A = B. Then r n — > ucp, so by Theorem 12.31 for all t > 0, 

X{t) = 1(0) + MB(t)) + k£ L 3 ,r - ^< r ) dW(s) = 0, 

where X(0) is the limit in probability of A n (0) — 6 n (0) as n — > oo, and $' r = V9 i r with 
$r(0) = 0. Thus, X(0) = 0, which shows that (i) holds. And as above, we obtain 
V^r = </?3,r = 0, which shows that (ii) holds. □ 

Let [A] denote the equivalence class of A under this relation, and let [S] denote the 
set of equivalence classes. If N = [A] G [S], then we define tpi^ = <fii,A, <^3,n = f3,Aj 
Xtv(0) = limA n (0), and 

T N (t) = I N (0) + ® N (B(t)) + kJ^ L 3 ,n - ^In) (B(s)) dW(s), (2.4) 

where &' N = tpi^ and $jv(0) = 0. Notice that by Theorem 12. 3^ if A r 1 , . . . , N m G [S] and 
AW G N k are arbitrary, then (B, A&, . . . , A^ m) ) (B,I Nl ,. . . ,2 Nm ) in D Rm +i[0, oo). 

It is easily verified that [S] is a vector space under the operations c[N] = [cN] and 
[M] + [N] = [M + N], and that N i— > X x is linear and injective. This gives us a one-to-one 
correspondence between [S] and processes of the form (12 .4j) . 

If X = f(B), where / G C°°(R), then we define N x = [A x ] G [S]. We may now identify 
X with N x , and will sometimes abuse notation by writing X G [S]. It may therefore 
be necessary to deduce from context whether X refers to the process f(B), the sequence 
A x = {A^}, or the equivalence class N x = [A x ]. Typically, there will be only one sensible 
interpretation, but when ambiguity is possible, we will be specific. 

Note that, using Q, we obtain <p ljX = /', ip 3X = & f", X x (0) = X(0) = f(0), and 
$x = / — /(0). Hence, by (I2.4p . we have T x {t) = X(t). Because of this, and because 
of the one-to-one correspondence between N G [S] and the process X/v(i) in (12. 4ft . we will 
sometimes abuse notation and write N(t) = N t = X x (t). Again, when there is a possible 
ambiguity as to whether N refers to an element of [S] or to the process In, we will be 
specific. 

2.3 The signed cubic variation 

If A G S, we define V A (t) = E}=M'(A„)) 3 and V A = {V n A }. Since 5 3 (V A ) = (S 3 (A n )f, 
it is easy to see from (12.21) that V A G S, fiyA = 0, and <£> 3 yA = (p\ A . Hence, if A = 6, 
then V A = V @ . We may therefore define the signed cubic variation of AT = [A] G [S] to be 
[V A ] G [S]. We denote the signed cubic variation of A" by [N]. We then have fiANj — 0, 
¥>3,[at] = <fil tN , and 2fljv](0) = 0, so that by ([22]), 

[Nj t =T m (t) = k f (<p 1>N (B(s))) 3 dW(s). 
Jo 

For example, suppose X = f(B), where / G C°°(R). Then [X] = {N x j. Since N x = [A x ], 
we have {N x j = [V aX ]. Note that V A * = {V AX } and 

[nt] [nt] 

V n AX (t) = $>,(A*)) 3 = Y^iXitj) - X^-r)) 3 . 

3=1 3=1 



s 



In other words, [X] is the equivalence class in S of the above sequence of sums of cubes of 
increments of X. By Theorem 12. 3[ {Xj t = Xpaif) is the stochastic process which is the limit 
in law of this sequence. Since <f% t x = f, we have f\\x\ — and <f3,ixj = (f) 3 , so that 

\X] t =l lxl (t) = K !\f\B{s))fdW{s). 
Jo 

In particular, taking f(x) = x gives {B} t = nW. 

2.4 The weak Stratonovich integral 

If A n , 9 n G S n , then we define 

\nt\ , . . k I \ 

(A n oe n )(t)^*^ W . 

If A, G S, then we define A o 6 = {A n o Q n }™ =1 . 

Lemma 2.5. If X = f(B), where f G C°°(R) and A G S, then A x o A G S. Moreover, if 
A = Q, then A x oAeA x o9. 

Proof. Clearly, A* o A n G S n and A x o A n (0) = for all n, so that Definition O^i) and (ii) 
hold. For a, b G K., we use the Taylor expansion 

/W t /(a) = f(x) + ±/" (x)(6 - a) 2 + i/ (4) (x)(fc - a) 4 + h(a, b)(b - a) 6 , 



4!2 

where x = (a + 6)/2, and |/i(a, 6)| < M(a,b) = sup xg [ aAb aVb ] |# (6) (z)|. For a derivation of this 
Taylor expansion, see the proof of Lemma 5.2 in [B]. 
Taking a = Bitj^i) and 6 = -B(tj) gives 

Ajffe_Q + A?(t 3 -) _ /(gfo-0) + /(gfo)) 
2 2 

= m) + ^m)A^ B + ^/ (4) (/3,)A4„ + 

where for each T, K > 0, there exists a finite constant Ct,k such that 

|-Ri,n|l{B*(r)<x> < Ct,k\ ABj tn \ 6 , 
whenever j/n < T. Choose ip?, G C°°(M) such that 

6j{A n ) = <p ltA (Pj)AB jtn + ^3,a(^-)A^„ + cp^AB^ + R j>n , 
where for each T, K > 0, there exists a finite constant Ct,k such that 

\Rj,n\l{B*(T)<K} < Cr,x| A_Bj jn | 7 , 
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whenever j/n < T. Then 

= (M, A )(^)Ai^ n + Q/'Vl,A + M,a) Wi) AB ln + K^ B ln + %,n, 

for an appropriately chosen smooth function h, and with Rj >n satisfying Definition I2.1( iii). 

It follows that A x o A £ S, and that fi^oA — ffi,A an d <^3,a x °a = f/'Vi.A + ff3,A- 
This implies that if A = 6, then A x o A = A x o 6. □ 

If X = /(£), where / £ C 00 ^), and N = [A] £ [5], we may now define XoN = [A x oA\. 
Note that if Y = g(B), where g £ C°°, and M E [S], then (A" + y)oJV = XoJV + yoJV 
and X o (AT + M) = X o N + X o M. From the proof of Lemma 12.51 we have 

1xon(0) = 0, (2.5) 

<Pl,XoN = ffl,N, (2.6) 
<^3,XoN = o/'Vl.iV + f<f3,N- (2.7) 



We may use these formulas, together with (12 .4p . to calculate T Xo n, given /, <pi t N, and <^ 3 , 



AT- 



We now adopt some more traditional notation. If X = f(B), where / £ C°°, and 
N £ [S], then 

J XdN = XoN E [S], 



and 

ct 



X(s) dN(s) = (XoN) t =Z XoN (t). 

As we noted earlier, there is a one-to-one correspondence between [S] and processes of the 
form (12.41) . We may therefore go back and forth between the above two objects according to 
what is more convenient at the time. We will use the shorthand notation dM = X dN to 
denote the equality M = J X dN. 

Before investigating our change-of-variable formulas, let us first consider some examples. 

Example 2.6. Let X = f(B) and Y = g(B), where f,g E C°°(R). Then 

J X dY = X o Y = X o N Y = [A x o A 5 
and A A ' o A Y = {A x oA Y }, where 

[nt\ 



(a? o aim = E " 2 ^ (a:) 



^ X((H) + X(t,) ro _ % ^ 
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In other words, J X dY is the equivalence class in S of the above sequence of Stratonovich- 
type symmetric Riemann sums. Also, J t X(s)dY(s) = TxoY(t), so that by Theorem 12.31 
JqX(s) dY(s) is the stochastic process which is the limit in law of this sequence. 

Example 2.7. Again let X = f(B) and Y = g(B), where f,ge C°°(R). Then 

J Xd[F]=Xo[F] = [A^o^ y ] 

and A x o V AY = {A x o }, where 



3=1 

= W xfe_ l) + x fe ) (yfe) _ yfe _ i))3 

In other words, J Xd[7] is the equivalence class in S of the above sequence of sums, and 
J *X(s)d[F] s = Xxo[Y\{t) is the the limit in law of this sequence. Recall that y?i,[y] = 
and <p 3> iYj = (g') 3 - Hence, by USD and (J277D, we have <pi t xo\Y\ = f<Pi,[Y\ = and 
<^3,Xo[y] = |/'Vi,[y] + f<P3,pr\ = fXg'Y, so that by (E3J, we have 

f X{s)d{Y\ s = k f f{B{s)){g\B{s))fdW{s). (2.8) 
Jo Jo 

Example 2.8. For one last example, let X = f(B), Y = g(B), and Z = h(B), where 
f,g,he C°°(K), and let iV = fY dZ. Then 

J X dN = X o N = X o [A y o A z ] = [A x o (A y o A z )], 

and 

(Aj o (A y o A z )„)(«) = (Af o (A^ o Af))(i) 

= ^(,- 1 + x fe )y fe -,H yfe ) (z( , |) . %i „, 

Hence, J X dN is the equivalence class in S of the above sequence of sums, and f*X(s)dN(s) 
is the limit in law of this sequence. 



3 Change-of-variable formulas 

We have already identified smooth functions of B with their corresponding sequences in S, 
as well as with their equivalence classes in [S]. In this section, it will be helpful to do the 
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same for J-^-measurable random variables, which can serve as initial values for the stochastic 
processes we are considering. 

Let T) be an J-^-measurable random variable, let A^(t) = r\ for all t > 0, and let 
A v = {A^}. Since <5j(A^) = for all j and n, we have that A v G S. We may therefore 
identify 77 with A G S, and also with N v = [A v ] G [S]. Note, then, that V2i >5? = </? 3)77 = 0, 
and rj(t) = jV"(t) = Z^t) = »7 for all t > 0. Note also that f X drj = 0. 

We begin with the following result, which tells us that every element of [S] has a unique 
decomposition into the sum of a smooth function of B and an integral against \B\ . 

Lemma 3.1. Each N G [S] can be written as N — T) + Y -\- V, where r\ is an J~^- 
measurable random variable, Y = g(B) for some g G C°°(R), and V = j 9(B) d{B\ for 
some 6 G C°°(E). 

Suppose N = rj + Y + V is another such representation, with Y = g(B) and V = 
J 0(B) d\B\ . Let c = g(0) — g(0) . Then rj=r] + c,g = g — c, and 9 = 9. In particular, there 
is a unique such representation with g(0) = 0. 

An explicit representation is given by rj = A r (0) = T N (Q), 9 = ip 3)N — ^<fi N and g chosen 
so that g' = <fi^ and g(0) = 0. 

Proof. Let N G [S]. Let 77 = N(0), 9 = tp^N — hfiN an d choose g so that g' = (fi^ and 
g(0) = 0. Let Y = g(B) and V = J 9(B) d[B]. To prove that N = r] + Y + V, it will suffice 
to show that 

N(t) = rj(t) + Y(t) + V(t) 

= N(0)+g(B(t))+ [ 9(B(s))d{B} s . 
Jo 

But this follows immediately from (" 12 .4[) and (12. 8p . 

Now suppose N(t) = ?j+g(B(t)) + f*9(B(s)) d{Bj s . Then E[N(t) \ Fg] = r) + g(B(t)) = 
rj +g(B(t)) a.s., which gives rj — rj + (g — g)(B(t)) = a.s. for all t > 0. Hence, there exists 
a constant c G M such that g — g = c, and it follows that rj = 77 + c. We then have 
M(t) = f*(6 - 9)(B(s)) dW(s) = a.s., so that E[M(t) 2 | JF*] = $1 \(9 ~ 9)(B(s))\ 2 ds = 
a.s. for all i > 0, which implies 9 = 9. □ 

We next verify that processes of the form V = j 9(B) d\B] behave as we would expect 
them to in regards to integration. 

Lemma 3.2. Let X = f(B), where f G C°°(R), and let 9 G C°°(R). If dV = 9(B) d{B], 
thenXdV = X9(B)d[B\. 

Proof. Let V = J 9(B) d{Bj, U = J X9(B) d{Bj, and N = J XdV. Since N(0) = U(0) = 
0, it will suffice to show that y2 1)( y = ipi^ N and (p 3j u = ^3^. By Example 12.71 tp\y = cpi t u = 0, 
f3,v — 9, and tp^u = f9. On the other hand, by (12. 6p and (12.71) . we have ipx,N — ffiy — 
and v?3,7v = |/'Vi,v + f<P3,v = f'9- □ 

We finally present our main result for doing calculations with the weak Stratonovich 
integral. 
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Theorem 3.3. Let N £ [S] and write N = rj + Y + V, where 77 is an Tf^-measurable random 
variable, Y = g(B), and V = J 9(B) d[B] for some g,6 £ C°°(M). Let X = f(B), where 
f £ C°°(R). Then 

J XdN = <5>(B) + ^ j{f"g'-f'g"){B)d{B\ + J X dV, (3.1) 
where $ £ C°°(R) is chosen so that $' = fg' and $(0) = 0. 

Remark 3.4. Since M = J X dN £ [S], Lemma 13.11 tells us that M has a unique 
decomposition into the sum of a smooth function of B and an integral against {B] . Theorem 
13.31 gives us a convenient formula for this decomposition. 

Remark 3.5. Theorem 13. 3} and the corollaries that are to follow, express equalities in the 
space [S]. Each side of (13.1 j) is an equivalence class of sequences of Riemann sums that 
converge in law. The equivalence relation is such that if we choose any sequence from the 
class on the left and any sequence from the class on the right, then their difference will 
converge to zero ucp. Note that this is a stronger statement than simply asserting that the 
two sequences have the same limiting law. 

Proof of Theorem E£B Since / X dN = f X dn + / X dY + / X dV and / X drj = 0, it 
follows from (12. 8p that we need only show 



X(s) dY(s) = $(B(t)) + - / (/V " f'g")(B(s)) dW(s). (3.2) 



By (j2~4"|) . we have 



t 

X(s)dY(s) 



o 



where M = X o Y. Recall that tp^y = g' and tp^y = 24 d'" ■ By fl2.6j) and (12. 7p . we have 
<Pi,m = fg' and (f 3>M = \f"g' + j^fg'"- Since $ M (0) = and $' M = = we have 
$az = $, and we also have 

1111 

ft Pit I Pill I £ t\H 

j// / I 1 j /// ^ /•// / ^ // L , w 

~ 8 24 ~ 24 61 ~~ 12 9 ~ 24 

= y 2 U"9'-f'g"), 

and this verifies (13.21) . □ 
Corollary 3.6. Let Y = g(B), where g £ C°°(R), and let <p £ C°° . Then 



ip{Y{t)) = <p(y(p)) + jf v/(y( a )) dY(a) - 1 jf ^'"(Y( S )) d[Y] s . (3.3) 
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Proof. Let X = <f'{Y) = f(B), where f — ip' o g. By Theorem 131 

JxdY = $(£) + 1 J (f"g' - /V0(B) d[B], 

where $ G C°°(R) is chosen so that $' = /gf' and $(0) = 0. Since {<p ° g)' = fg', we have 
$ = (<pog)- (<fog)(0). Also, 

fg' - fg" = ((</' o g)(g') 2 + W o g)g")g' - (</ o «?)«/</' = (</' o <?)(</) 3 . 

Thus, 

-t rt 



[ <p'(Y(s))dY(s) = f X(s)dY(s) 
Jo Jo 



1 " 



(<p o g)(B(t)) - (<p o g)(0) + - I &"og)(B( S ))(g'(B( S ))fdlB] 



cp(Y(t)) - <p<y(p)) + ip"'{Y{s)){g'{B{s))? dW(s) 



By ([22]), this gives 

jT ^CrW) dy(«) = ^(y(t)) - ^(y(o)) + ^ jT d[Y] s , 

which is (Q. □ 

Corollary 3.7. Lei iV G [5] and write N = rj + Y + V , where r\ is an F^-measurable 
random variable, Y = g(B), and V = J 9(B) d[B] for some g, 9 G C°°(R). Let X = f(B) 
and Z = h(B), where f,he C°°(R). Then: 

If dM = X dN, then Z dM = ZXdN - ~{fg'ti){B) d{B\ . (3.4) 

Moreover, the above correction term is a "weak triple covariation" in the following sense: If 
V = {V n }, where 

\nt\ 

i=i 

then V eS and [V] = f(f'g'h')(B) d{Bj. 

Proof. Let N, X, and Z be as in the hypotheses, and let M — J X dN. By Theorem 13.31 

M = + 1 J (fg' - fg"){B) d\B\ + J X dV, 

where $ G C°°(R) is chosen so that $' = fg' and $(0) = 0. Hence, by Lemma 

J ZdM = J Zd$(B) + iy (fg'h - f'g"h){B)d{B\ + J ZXdV. (3.5) 
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J Zdf(B) = ?(5) + i J(h"&-hl&)(B)&\Bl 
where * G C°°(E) is chosen so that = /i<3>' and $(0) = 0. Theorem E3] also gives 
J ZXdY = ^(B) + ± f{{fh)"g'-{fh)'g"){B)d{Bl 

where ^ G C°°(IR) is chosen so that \&' = fhg' and \&(0) = 0. Note, however, that this 
implies ^ = ty, which gives 



Zd<S>(B) = JzXdY + ^ j{h"&-h'&'-{fh)"g'+{fh)'g"){B)d{Bl 
Substituting $' = fg' into the above and simplifying gives 

J Zd<$>{B) = JzXdY + ^J (f'g"h - f'g'h - 3f'g'h')(B) d[B]. 
Substituting this into (13. 5p gives 

J ZdM = J ZXdY-^ J {fg'h'){B)d{B\ + J ZXdV 
ZXdN-\ J{fg'h'){B)d{Bl 



and this verifies ( 13.41) . 

Finally, if V = {V n }, then 5j{V n ) = 6j(A% )^(A^)^(Af ). From (Q, we see that V G 5, 
</>i v = 0, and 993 y = <£>i x<^i v<^i z = f'g'h'. Since [V]o = 0, it follows from Example 12 . 71 that 
[V]' = /(/W)(5)d[fi]. □ 
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